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Abstract. The degree of entanglement of random pure states in bipartite 
quantum systems can be estimated from the distribution of the extreme Schmidt 
eigenvalues. For a bipartition of size M > N, these are distributed according to a 
Wishart-Laguerre ensemble (WL) of random matrices of size X M, with a fixed- 
trace constraint. We first compute the distribution and moments of the smallest 
eigenvalue in the fixed trace orthogonal WL ensemble for arbitrary M > N. Our 
method is based on a Laplace inversion of the recursive results for the corresponding 
orthogonal WL ensemble by Edelman. Explicit examples are given for fixed N and 
M, generalizing and simplifying earlier results. In the microscopic large-N limit 
with M — N fixed, the orthogonal and unitary WL distributions exhibit universality 
after a suitable rescaling and are therefore independent of the constraint. We 
prove that very recent results given in terms of hypergeometric functions of matrix 
argument are equivalent to more explicit expressions in terms of a Pfaffian or 
determinant of Bessel functions. While the latter were mostly known from the 
random matrix literature on the QCD Dirac operator spectrum, we also derive 
some new results in the orthogonal symmetry class. 
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1. Introduction 

The theory of random matrices finds appfications to the most diverse physical 
situations, and Wigner and Dyson are usually referred to as the pioneers of this field 
for their works on nuclear spectra. However, many years before their seminal papers, 
John Wishart had already introduced random covariance matrices in his studies on 
multivariate populations [1]. The Wishart ensemble (also named Laguerre or chiral, 
and called WL hereafter) contains random N x N covariance matrices of the form 
W = X^^X where X is a M x matrix with i.i.d. Gaussian entries (real, complex or 
quaternions labelled by the Dyson index f3 = 1,2,4 respectively). The joint probability 



density function (jpdf) of its TV non-negative eigenvalues {A^} is given in Eq. (2.8). 
A related ensemble of random matrices is the fixed-trace Wishart-Laguerre (FTWL), 
which contains WL matrices with a prescribed trace TrCW) ~ t > and whose 



eigenvalue jpdf is given in (2.10). Recent results about the density of eigenvalues 
of FTWL matrices can be found in [5J [31 The study of fixed and bounded trace 
ensembles of non-chiral random matrices has however a longer history [5] . In [51 [7] the 
spectral density of the related, fixed trace Gaussian Unitary Ensemble was computed 
explicitly for finite N, while for works on non-chiral fixed trace /3-ensembles we refer 
to [El [9]. 

Two recent applications of WL and FTWL ensembles that motivate our study are 
effective theories of Quantum Chromodynamics (QCD) and QCD-like theories, where 
we refer to [TO] for a most recent review and references, and the statistical theory of 
entangled random pure states in bipartite systems, see [IT] for an excellent review, 
respectively. 

Entanglement is indeed one of the most distinctive features of quantum systems, 
and it is a crucial resource in quantum computation issues, as performances of 
quantum computers will heavily rely on the possibility of producing states with 
large entanglement [12 US]. In order to quantify the degree of entanglement of a 
given quantum state, it is useful to introduce an unbiased benchmark of states with 
the lowest degree of built-in information: pure bipartite systems (defined below) 
with a Hamiltonian of size M ■ N constitute a typical example where well-behaved 
entanglement quantifiers can be defined, such as the von-Neumann or Renyi entropies 
of either subsystem (13] , the so called concurrence for two-qubit systems ^14j or other 
entanglement monotones [TS] [16] . 

The main focus of this paper is on the cumulative distribution and density of 
the smallest FTWL eigenvalue and its relation to the corresponding quantity for the 
unconstrained WL ensemble: these two distributions arise naturally in the two settings 
described above (entanglement and QCD, respectively). The upshot is that the smallest 
Schmidt eigenvalue (a relevant quantity in the entanglement setting discussed below) 
is precisely given by the smallest eigenvalue of one of the FTWL ensembles with 
Tr(W) = t = 1. Which of these ensembles to choose depends on the global symmetry of 
the quantum system, and we will report on the orthogonal (/3 = 1) and unitary (/3 = 2) 
symmetry classes, but not on the symplectic {/3 ~ 4) class. Another application of 
the smallest FTWL eigenvalue is related to the Demmel condition number ,17 , which 
quantifies how hard certain numerical problems about matrix inversion are. 

The results of this paper can be grouped into two parts: 

(i) First, we study the full distribution of the smallest eigenvalue of FTWL orthogonal 
ensemble {f3 = 1) at finite N, M . 

(ii) Next, we take the large- TV, M limit of the smallest eigenvalue with v ~ M — N >Q 
fixed (the so-called microscopic limit) for both /3 = 1,2. 

Let us summarize first which results were previously known for both FTWL and 
unconstrained WL: 
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• For FTWL ensemble at finite N, M, the full distribution of the smallest eigenvalue 
was derived for /3 = 1, 2 and M = iV in [18jU and for /? = 4 in [IT]. In [15], 
a conjecture by Znidaric [22) was proven for the first moment. For general real 
/3 > and m e N, where m = (/3/2)(A/ - + 1) - 1, Chen et al. gSj reported 
formal expressions for the smallest eigenvalue cumulative distribution and density 
in terms of cumbersome sums of Jack polynomials. More explicit expressions are 
given for the case /? = 2 in terms of finite sums. 

• For WL ensemble at finite N, M, the first results for the smallest eigenvalues of 
the orthogonal WL ensemble go back to Edelnian JH [25] who gave a recursive 
prescription. Curiously, analogous recurrence relations for /3 = 2,4 have not 
been worked out to date, while for /3 — 2 a, closed expression exists in terms 
of determinants [IS]- These results were extended by Forrester [57] to the case of 
real /3 > for special values of = M—N. Forrester again waived some limitations 
on 1/ and further generalized these results [28]. The most up-to-date and general 
formula for the smallest eigenvalue distribution involves hypergeometric function 
of matrix argument (HFMA), see eq. 2.13b in [55], as well as our Appendix |Aj 
Note that for (3 = I this requires v to be odd. For /3 = 1,4 expressions in terms 
of Pfaffians were reported [2S^, with the same restriction for (3 — 1, and for /3 = 4 
the values of i/ > are half- integers. 

The next natural step is to take the large N, M limit for the smallest eigenvalue. For 
real f3 > keeping v fixed, this was done in [35] for the WL and subsequently in [33] for 
FTWL ensembles, expressing them in terms of HFMA. It turned out that in this large 
TV limit (after a suitable rescaling) the fixed trace constraint does not matter, and in 
that sense the results are universal (for a different limit, namely c = N/M fixed, see 
the recent paper [50]). 

The issue of universality in the fixed (or bounded) trace ensembles had been 
addressed earlier in the literature, mainly for the non-chiral ensembles. It was shown 
in [6] that in the macroscopic large N limit where the oscillations of the density are 
smoothed, the semi-circle or its generalizations can be matched in the constrained 
and unconstrained ensembles. In contrast, the higher order connected correlators 
become non-universal when considering non-Gaussian generalized WL ensembles, with 
a polynomial potential in W in the exponent. On the other hand, in the microscopic 
limit in the bulk of the spectrum, where the oscillatory behavior of the density is zoomed 
into, the constraint was found to be immaterial and the universality of the sine-kernel 
was established in |31j. The same result was provided later in a mathematically more 
rigorous way for non- invariant generalizations of WL in [331 [33] . 

Because of this universality (i.e. constraint-independence), it is extremely useful to 
recall the large N results derived independently for the smallest eigenvalue distributions 
in the application of unconstrained WL ensembles to QCD [SH [35l [36]. Here also 
non-Gaussian generalisations of WL were considered, and determinantal or Pfaffian 
expressions of Bessel functions were derived for /3 = 2 [55] and /3 = 1 [35] , respectively. 
In fact much more general results were derived there, including an arbitrary number 
of characteristic polynomials of random matrices (so-called mass terms) in the weight 
function. In the second part of this paper we will show how these results in the QCD 
literature and the aforementioned ones in terms of HFMA are related. 

The presentation of the paper is organized as follows. 
In order to be self-contained in section [2] we provide some background material about 



the relation between bipartite entanglement and FTWL (subsection 2.1 ). In subsection 



J In studying the distribution of the Demmel condition number, Edelman 1191 has derived the same 
formulae in a completely different setting. An exact mapping between the two problems is possible, 
but nontrivial 12011211. 
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|2.2| we give definitions and notations used for WL and FTWL ensembles. The reader 
familiar to either or both of these topics may skip the corresponding subsection(s) . 

In the next section [3] we derive our new results for FTWL ensembles at /? = 1 for 
arbitrary finite N and M, both for odd and even v. There we briefiy recall the results 
for standard WL by Edelman on which we heavily rely. This section also contains our 



results for arbitrary moments in subsection |3.3| and numerical checks in 3.4 



Section |4] brings us to the universal microscopic limit for fixed v = M — N . It 
provides equivalence proofs between the known results from the QCD Dirac spectrum 
literature, given in subsection |4.2[ and the very recent results in terms of HFMA, given 



in subsection 4.1 In subsection |4.3| we explicitly compute the large N limit of the scaled 
smallest eigenvalue for (3 — 1 and u = 0,2 from their finite N expressions and confirm 
the universality results. For (3 — 2 the above mentioned equivalence is established in 
subsection |4.4[ and for (3 — 1 in subsection [43] including new results for the cumulative 
distribution and for even v = 2. We also report on the general scaling for moments in 
subsection 4.6 before concluding. 

Some technical details about the definition of HFMA and a universality check for 
(3 = 1 and = are deferred to the appendices. 



2. Bipartite entanglement and FTWL ensembles 

2.1. Bipartite entanglement 

Let 5f(^*^) be a TV • M-dimensional Hilbert space which is bipartite as 5f(^^^) = 
:k'x^ ^ J£^*^\ where N < M without loss of generality. For example, A may be 
taken as a given subsystem (say a set of spins) and B may represent the environment 
(e.g., a heat bath). Let {|«'^)} and {|q!^)} be two complete bases of ^[^^ and 
respectively. Then, any quantum state \ip) of the composite system can be decomposed 
as: 

N M 
i—l a — 1 

where the coefficients Xi^a's form the entries of a rectangular (A'^ x M) matrix X. 

In the following, we will consider entangled random pure states \ip). This means 
that: 

(i) I"!/)) cannot be expressed as a direct product of two states belonging to the two 
subsystems A and B. 

(ii) the expansion coefficients 

X 2. 3)1" G random variables drawn from a certain 
probability distribution (see below). 

(iii) the density matrix of the composite system is simply given by p = IV') ('/'I with the 
constraint Tt[p\ = 1, or cquivalently {^(14^) = 1. 

We will not consider statistically mixed states here, and we refer to [37] and 
references therein for recent results on mixed states. 

The density matrix p of a quantum state is a very important quantity as it allows 
to compute expectation values of observables. For an entangled pure state \ip) of a 
bipartite quantum system it can then be straightforwardly written as: 

^ = E E 4^ 1*^) I ® I"'') (^"^ I ' (2-2) 

where the Roman indices i and j run from 1 to and the Greek indices a and /3 run 
from 1 to M. 
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In some applications, it is useful to separate the contribution of the subsystem A 
under consideration from the environment B. Expectation values of observables O^i on 
the system A can be obtained by "tracing out" the environmental degrees of freedom 
(i.e., those of subsystem B) and defining the reduced density matrix pA = TtbIp] as: 



M 

PA = TTB[p]^y(a^\p\a^). (2.3) 



Using the expansion in Eq. (|2.2[) one gets 



N M N 

= E E^^"4" - E (2-4) 

where Wij 's are the entries of the N x N covariance matrix W = XX^ . 

Proceeding further, one could also obtain the reduced density matrix pB — Tr^[p] 
of the subsystem B in terms of the M x M matrix W = X^X and find that W and 
W share the same set of nonzero (positive) real eigenvalues {Ai, A2, . . . , Aat}, called 
Schmidt eigenvalues. 

In the basis of eigenvectors of W, one can express pA as 

N 

PA = J2K\Xt){Xf\ (2.5) 

1=1 

where |Af )'s are the normalized eigenvectors of W = XX^^. The original composite 
state in this diagonal basis reads: 

N 

H)=J2^/)H\Xt)^\Xf). (2.6) 



Eq. (2.6) is known as the Schmidt decomposition, and the normalization condition 



{ip\ip) = Ij or cquivalently Tr[(0^] = 1, imposes a constraint on the eigenvalues, 

In the Schmidt decomposition (2.6), each state |A^) |Af ) is separable, but their 
linear combination \ (depending on the set of Schmidt eigenvalues) cannot, in general, 
be written as a direct product IV') = \4>^) \(t>^) of two states of the respective 
subsystems, i.e. it is entangled. Knowledge of the Schmidt eigenvalues {Ai, A2, . . . , Ajv} 
of the matrix W is therefore essential in providing information about how entangled a 
pure state is. 



For random pure states, the expansion coefficients in eq. (2.1 1 can be typically 
drawn from an unbiased (so called Hilhert- Schmidt) distribution 

Prob[X] cx 5 (Tr(XX^) - l) . (2.7) 



The meaning of eq. ( |2.7[ ) is clear: all normalized density matrices are sampled with 
equal probability, which corresponds to having minimal a priori information about 
the quantum state under consideration. This in turn induces nontrivial correlations 
among the Schmidt eigenvalues (which are now real random variables between and 
1 whose sum is 1) and makes the investigation of several statistical quantities about 
such states quite interesting. The jpdf of Schmidt eigenvalues for a Hilbert-Schmidt 
distribution of coefficients was derived in 38 and turns out to be exactly of the FTWL 



form (2.10) with t = 1, where the Dyson index /? = 1,2 corresponds respectively to 
real and complex X matrice^ The delta constraint there indeed guarantees a proper 
normalization of the reduced density matrix Tr[p^] = 1. 



§ These two cases in turn correspond to quantum systems whose Hamiltonians preserve (/9 = 1) or 
break (/3 = 2) time-reversal symmetry. 
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Why is the smaUest Schmidt eigenvalue distribution interesting at aU? We first note 
that due to the constraint J2iLi -^1 = 1 and the fact that all eigenvalues are nonnegative, 
it follows that 1/N < A,nax < 1 and < Amin < 1/N. Now consider the following 
limiting situations. Suppose that the smallest eigenvalue Amin — niini{Ai} takes its 
maximum allowed value Then it follows immediately that all the remaining 

{N — 1) eigenvalues must be also identically equal to Xi = 1/N. In this situation eq. 



(2.6) tells us that \ip) is maximally entangled. On the other hand, if Amin = (i.e., it 
takes its lowest allowed value) or close to 0, while this will not provide much information 
about the degree of entanglement of ji/'), it actually tells us that one component in the 



Schmidt decomposition (2.6) can be safely ignored. In other words, the 'effective' 
dimension of the Hilbert space !Ha has been reduced from N to N — I. The proximity 
of the smallest eigenvalue to zero, therefore, gives information about the efficiency of 
this dimensional reduction process. 

For more references on entangled random pure states we refer to [llj . 

2.2. WL Ensembles with and without fixed trace 



The jpdf of non-negative eigenvalues of the unconstrained WL ensemble is given by: 



T(wi^)(Ai,...,A^) = <i,exp 



2 

1=1 



|(1 + M-7V)-1, . 



1=1 



where J^^^f is a known normalization constant: 



/ a\ |JVA/ N r ( 1 + ^ 



..r(i + fj)r(f(M-« + i)) 



and A(A) = n^xfel-^j ~ ^k) is the Vandermonde determinant. 

On the other hand, the jpdf of the eigenvalues A; e [0,t] of the FTWL ensemble is 
given by: 

T(^^)(Ai, ...,X^;t) = &^\,{t)5 E - n A,?^^^-^+^^- V(A)r . (2.10) 

\i=l / i=l 

For i = 1 it coincides with the distribution of Schmidt eigenvalues, where the delta 
function guarantees that Tr[p] = 'Yli^=i — ^- The normalization constant for t = 1 
reads [ID]: 

r(« -r^fi) ni- r(Mjv/3/2)(r(i + ;3/2))^ 

The presence of a fixed-trace constraint is known to have crucial consequences on 
(connected) spectral correlation functions, both for finite- iV and in the macroscopic 
large- iV limit fBlIM]- However, in the microscopic large- TV limit the correlations become 
universal, and we shall exploit this in the next section. Let us first introduce the crucial 
quantities for this paper, the cumulative distributions qN,u (x) = Prob[Amin > x] (also 
called gap-probabilities) , and the corresponding densities for the smallest eigenvalue of 
both ensembles: 

N / N \ ^ B 

<i^Z\--^t) - C%it) / \[dX,5[Y^X.-t]\[xf''^''''\l:.{X)\f^ , (2.12) 
J[xM" t=i \7~{ ) 7=1 



pK?(^;0 = -^C^(^;t) , (2.13) 
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N N 



.r'(-) =<m/ n^^^^-^^-^''n^f'"'^">wr (2-14) 

pZ\-) (2-15) 

where the superscript distinguishes fixed-trace (FT) and ordinary WL. We have the 
following normalization conditions: 

OO nt/N 

dx p^^^^Hx) = y dx p^^l\x; i) - 1 , (2.16) 

iT(0) = 'z!vT(0;i) = i. (2.17) 

The tt\v moments of the smallest Schmidt and WL eigenvalue are then given by 

(AL„)^7,^ / dx x'p'^^l\x) , (2.18) 

/•OO 

(ALJK- / dxx'p^Z^\x). (2.19) 







where we define p)^ J{x) p^q J{x; 1) for later convenience. All density correlation 
functions and thus also the smallest eigenvalue distribution in ensembles with and 
without fixed trace can be related via an inverse Laplace transform, see e.g. ^Slj [ji] For 
the smallest eigenvalue distribution this relation reads: 



£[p!v7(^;i)](-^)= / dtp'-^l'{x-t)e 







(^) pivT^(2..) . (2.20) 



'N.M 



This is the main technical tool we use in the next section. For finite N, M we are 
going to use known explicit expressions for WL smallest eigenvalue densities in order 



to derive the sought density for FTWL via the inverse Laplace transform of eq. (2.20) 



3. New results for FTWL ensemble at finite N, M and (3 — 1 

In this section we will focus on the (3 = 1 ensemble only. Our strategy is very simple. 
We will take the known results for the standard (3 = 1 WL ensemble at any finite N, M 



derived by Edelman [25] and invert the Laplace transform in eq. (2.20). The same 
strategy was followed by Chen et al. [23] . however, they applied the inverse Laplace 
transform directly to the HFMA result by Forrester [21] valid for real /? > 0. This 
led to expressions that are formally exact, but not very transparent or manageable. 
Furthermore, in contrast to [531 HH] we will not be restricted to odd = M — N ioi 
(3 = 1- The special case for i/ = was derived independently in [TS] and [15] , 

3.1. Odd V = M - N 

We begin with the simpler case. Edelman [25] states that the smallest eigenvalue 
distribution for a WL ensemble with (3 = 1 and v = M — N odd can be written as 
follows 15] 



Pn 



(WL) 



'^\x) = 2f -lc^.,.:r(^-l)/2e-^-/2V.(x) (3.1) 



II For this reason we have kept the fixed trace to be a free parameter t in the jpdf ( 2.10 i. 

% Compared to the notation in 25 we spell out the parameter p explicitly contained in cjv v there. 
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where h^^^^x) is a polynomial of degree {N — l){h' — l)/2 and cn,!^ is the following 
constant: 



^2-JV^/2r((Ar + i)/2) J\ T{j/2) 
CN.. -j= 11 



^^T{{N + j)/2) 



As such, hi^,y(x) can be written as: 

(JV-1)(^-1) 



(3.2) 



(3.3) 



with some rational coefficients q^, which depend also on N, v. The polynomial hN,i/ is 
determined via a simple recursion relation explicitly given in . In addition, a simple 
Mathematica® code is given in Appendix A of f25| to generate the hN^^{x). The first 
two polynomials read: 

and so on using the given recursion. In hA, l'C\x) is a generalized Laguerre 
polynomial 

Now, it is easy to observe that (no matter what the coefficients qfc are) the form 



(3.1) lends itself to a very friendly Laplace inversion. More precisely, take the inverse 
Laplace transform of the fundamental relation ( 2.20 1 after inserting eq. ( 3.3 ) into (3.1). 
One obtains: 

(JV-l)(t.-l) 



C 



(1) 



N.M 



X 



k=0 



^2s-^(2k+u+l~N{N+u))/2^-Nsx ] 



(3.5) 



Computing the inverse Laplace transform and setting t = 1, one gets the final general 
formula: 

(FT), , q, 2^2k+.+l-NiN+.))/2 

^^.''^^^-kW/ r((Ar(7V + .)-2fc-.-l)/2)^ 

X (1 - „ ^^)_ (3 6) 

where 6{z) is the Heaviside step function. This equation is our first main result of this 
section. The algorithm to compute pj^P (2:) for odd v works as follows: 
(i) Compute the polynomial hiq_y{x) for the sought (A^, v) using Edelman's recursion 
relation [23; 



(ii) Extract the coefficients of the polynomial from (3.3); 

(iii) Insert these coefficients in the general formula (3.6 1. 

To give some explicit examples we have worked out in all detail the cases v = 1,3, 
which are based on eq. (|3.4|): 

(FT) / N 



A^C^'^^+i2-^(^+i)/2(l - Nx)-^+^^^+^ 



)/2 



K 



(1) 

N.N- 



^r(-l + A(A + l)/2) 



(1 - A^a;) 



(3.7) 



(FT) / N 



r(3 + A^) 
^(1) 2(A + 1)(A/-- 1)! 

22+k-N(N+3)/2 

^ T{-2-k + N{N + i)/2) 



k=0 



r(4 + fc)fc! 



a;i+'=(l- A^x)-3-'^-+^(^+3)/2 



(3.8) 



Note a typo in Lemma 4.5. in [25| in the representation of Li^\—x) 
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We stress, however, that the general formula (3.6) and the algorithm above provide 



an explicit and user-friendly solution to the problem of computing Pfq J{x) for any 
desired value of {N^v). The obtained results are much more explicit and manageable 
than those expressed in terms of a finite sum over partitions in [23] . To illustrate our 
algorithm we have generated examples with higher v. 

We wrote a simple Mathematica® code that generates explicit expressions for the 
smallest eigenvalue density (the code is freely available on demand) . A plot for different 
values of odd v and = 9 is provided in fig. [T] 



400 




Figure 1. Smallest eigenvalue density p^^' {x) for N = 9 and u = 1, 3, 5, . . . , 13 
odd, using formula \3.6\ and a simple Mathematica® implementation of our 
algorithm above. The running time for the full evaluation was less than 10 sec. on 
a standard laptop. 



3.2. Even v = M ~ N 

Now, we turn to the more complicated case v = M — N >Q even. Also in this case, we 
can derive a general formula for p]^ ^ {x) involving the same polynomial coefficients as 
in Edelman's recursion. To the best of our knowledge, apart from u = Q |18| . no such 
expressions were previously known for the FTWL ensemble. 

Edelman [53] states that the smallest eigenvalue distribution for a WL ensemble with 
(3 — 1 and v even can be written as follows: 

p^;^^\x) - 2-hN,.x^''-'^/^e-''-/^[fNA^)Uiix) + gNA^)U2ix)] , (3.9) 

where fN,v{x) and gN,ij{x) are polynomials of degree at most (A^ — l)v /2, and cjq.u is 
given in eq. (3.2). The two related functions Hi (a;) and 1X2(2;) are defined as follows: 

'A^- 1 1 



lli(x) = U 



2 ' 2' 2 



ll2(x)^-^u(^^,-,-j= Ui(x)', (3.10) 

where ?7(a, h, z) is a Tricomi confluent hypergeomctric function given by: 

1 1"°° 

U(a, b,z)^--- dre-^W^-^T + l)-"+''-i . (3.11) 

r(a) Jo 
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Let us now write the polynomials fN,u{x) and gN,v{x) as: 
f(JV-i) 
fnA^) = X! ^^^^ ' 

f(JV-l) 
fe=0 



(3.12) 



(3.13) 



with some rational coefficients ffe and Qk^ which depend also on v. As in the odd v 
case they follow from a recurrence relation given in [25^ along with a Mathematica® 
code to generate them. The first two examples read: 

r(iv)y^ 



fN.u=o{x) 



2A^-ir(f)r(^) 



gN,u=oix) = 



f (^\- ^ 2 j M 2 j r 2 , . 



Applying now eq. (2.20) to (3.9), we get: 

f(JV-l) 



c 



(1) 



(3.14) 



(3.15) 



(3.16) 



N.M 



where 



*,(^,s) ^ x^+(-i)/^e--^"(2g)^-+i+ '^'""'"^'^' Ui(2ga;) 

rfr T^(r+ l)~f-i(2s) 







*,(x,s) ^ ^ ^fc+(.-i)/2^-.iv.(2,)fc+i+ -^-^'"+-' 1x2(2.0:) 



4 r I' ^+i ' 



dT r^(T + l)-T-i(2s)'^+ 



^ + l-iV(W + ^) 



-san(Af+T) 



Here we have inserted the integral representation (3.11). 

Taking the inverse Laplace transform of the quantities in square brackets, and setting 
t = \ we can eventually write: 



^(1) f(JV-l) ^ + l + 2fc-JV(W + ,.) 

pSPW = ^2-^c^,. ^ -^^^^^^^^-^[f,0,(a:)-0,H,(x)] , (3.17) 



=n r 



fc=o 



where: 

0fc(x) 



a; 



fc + (i^-l)/2 /-oo 



r(^) 7o (r + l) 



W-3 

T 2 iVCAr+M)_M_3_2fc 

dr ^(1 - x{N + r)) 5 0(1 _ a;(7V + r)) 



JV(jV+i/)-i^-l-2A; 
2 



( 



jV(iV+i/)+jV-t^-2-2fc 



2;'=+— (1 _ f^x) 



N(N + iy)-iy + N -4-2k 



f N-l N -2k + {N -1){N + iy + 2) l\ 



X 2F1 

4 r{^: 



(3.18) 



dr T 2 (r + 1)" 



2 - X 
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JV(JV + i/)-i^-3-2fc 

X (1 - x{N + r)) 5 ^(l - x{N + t)) 



N - 1 



jV(jV+t/)-iv-l-2fc 
2 



^(1 - Nx} 



~k+ 



(JV-l)(JV + ^ + 2) 



X 2F1 



N(N+v)+N-v-2k 
2 

N +1 N -2k -iy + N(l + N + iy) 1 

, hi: ^ --.N 

2 ' 2 ' 2 ' X 



(3.19) 



In the second step we have performed the integrals and apphed some simple algebra, 

and 2-^1(1, b] c; z) is a standard hypergeometric function. 

Eq. (3.17) together with eqs. (3.18) and (3.19) is the second result of this section. 

/t"l' 

Again, the algorithm to compute p\s[ y foi' even v works as follows: 

(i) Compute the polynomials fN,v{x) and gN,iy{x) for the sought {N,v) using 
Edelman's recursion E51: 



(ii) Extract the coefficients % and of the polynomial from (3.12) and (3.13); 

(iii) Insert these coefficients in the general formula (3.17). 



To provide explicit examples we give the expressions for = 0, 2 based on eqs. (3.14 ) 

eviously derived in |18| : 

^-Ar/2(^_^^)(iV^+Ar-4)/2 ^ 



and (3.15), where the result for v = Q was previously derived in 

ivr(iv)r(ivV2) 



(FT) I 



2W-ir(iV/2)r((7V2 + 7V_ 2)/2)' 

;N + 2 N -1 N^+N -2 

X 2F1 



N - 



1 



and 

where: 

0^(2; 



C^N^N+2 TUN + l)/2) 



K 



(1) 

N,N+2 



27r 



(3.20) 



(3.21) 



N-l 

2(3-JV(Af+2))/2^(l-JV)/2^^ _ ^^^-3+3^/2+^^2 ^ 

fc=0 

3A^ 



2'=(iV + 1)! 



X 



X 



1- Nx 



N-l N N"^ 



fc!(iV-l-fc)!(fc + 2)! 
1 



-2-fc:iV- 
2 ' X 



(3.22) 



N-2 



i>j,{x) = 2(^-NiN+2))/2^il-N)/2^^ _ j^^y3+3N/2+N^/2 ^ 



2'''(iV + 1)! 



1- Nx 



N + 1 



N ^ N^ 



k=0 

3N 



fc!(iV-2-fc)!(fc + 3)! 



k:N ■ 



(3.23) 



2 2 2 2 X ^ 

In order to illustrate our algorithm above we again wrote a simple Mathematica® 
code that generates explicit expressions for the smallest eigenvalue density. A plot for 
different even > and iV = 9 as an example is provided in fig. [2] 

3.3. Moments for odd v 

Arbitrary moments (A^j;,j)^P^^-^ can be computed easily in a closed form for odd v 
from eq. (3.6), using the following integral formula: 



dx 



\\-Nxf'^ =N-'^V,{uj,i) 



(3.24) 
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0.0] Offi a.03 OM 



Figure 2. Smallest eigenvalue density P^^^ (a;) for N = 9 and i/ = 2, 4, 6, . . . , 12 
even using formula and a simple Mathematica® implementation of the 

algorithm above. The running time for the full evaluation was less than 50 sec. on 
a standard laptop. 



Here B(ix),^) is Euler's Beta function. We thus obtain 



(A 



I \(FT) 

min/Ar,^=2m+l,^=l 



c 



K 



(1) (Af-l)m 9m+l + (2fe-Af(7V+y))/2 

^ ^ CN,2rn+l 2^ 



(1) 



N,M 

xB(£ + fc + l + TO, 



fc-TO-l|. (3.25) 



fc=0 ^( 2 

7V(7V + 2m + l) 
2 

To give an example, the £th moment for v = 1 is given as follows 
In particular for the first moment [£ = 1) or average value we get the following answer 

(FT) _ 2 



(X^ \(FT) ^ r,~N(N+l)/2 ^(1 + tj 

K min> N,^=l,P=l ^ mr{e + N{N +l)/2) 



{Xn 



For comparison at = the large- behavior is as follows 

{X.nin}^rJ=0,P=l - ^/N' for TV » 1 , 

where: 



for > 1 



c = 2 



1 - 



-erfc(l 



0.688641 



(3.26) 

(3.27) 
(3.28) 



(compare with subsection 4.6). Although the computation of moments for even v is 



possible, based on our explicit expressions given previously, we did not find short closed 



expressions as in eq. (3.25) above 



3.4- Numerical checks 



In fig. [3] we compare a few theoretical densities from the previous subsections 3.1 and 
3.2 for finite N with the corresponding numerical results. These are obtained as follows 
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(i) we generate n « 10*^ real Gaussian M x TV matrices X (where N — 7 ,M — N + v). 

(ii) for each instance we construct the Wishart matrix W — X^X. 

(iii) we diagonahze W and collect its N real and non-negative eigenvalues {Ai, . . . , Xn}- 

(iv) we define a new variable < ^.i < 1 as ^ii = Amin/X^i^i 

(v) we construct a normalized histogram of /ii. 

The agreement between theory and simulations is excellent. 




v=1 NUM 
v=2NUM 
v=3 NUM 
v=1 THEO 
v=2 THEO 
v=3 THEO 



Figure 3. The density p^^'^ ^{x) of the smallest Schmidt eigenvalue from 
numerical simulations compared with our predictions eqs. and 
for = 1,2,3 respectively (solid lines). 



4. Equivalence proofs and nevif results in the large- A^, M limit 

We now turn to the microscopic large- A^, M limit for the orthogonal and unitary 
ensembles (/3 = 1,2), with v = M ~ N > Q fixed. More precisely, let us define 
the following microscopic scaling limits at the origin: 

Notice that the mean level spacing at the origin (the TV-dependent rescaling factor) 
is different for FTWL and standard WL ensembles. This fact was also observed when 
comparing the corresponding macroscopic densities in the Gaussian ensembles with 
and without constraint [B]. 

It has been shown that in the microscopic limit (i) the smallest eigenvalue 
distribution in non-Gaussian generalizations of WL ensembles is universal |351 136] . 
and (ii) for all real /3 > 0, the microscopic limit for the smallest eigenvalue distribution 
is the same for WL and FTWL ensembles in all cases where a representation in terms 
of HFMA exists [23]. We also recall that universality (ii), that is the independence of 
the constraint, was shown earlier for all microscopic density correlation functions in 
the bulk of the spectrum for unitary non-chiral ensembles for monic even potentials 
[5T| or i.i.d. matrix elements [5^ [55] , 
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Building on [23], from now on we shall assume that: 

Q'J^'^^y) = QrHy) (4.3) 

p(WL)(^)^p(FT)(y)^ (4.4) 

for all /3 and y. Therefore, there is no reason to keep the superscripts and 

We will thus use the unified notation Q^f (y) and Pl'^\y), labeled by /? onljQ Note 

that one still has 

Pl^\y) = -Qi^\y)' , (4.5) 

dyPi^\y)^l, (4.6) 



Q(«(0) = 1 . (4.7) 

What are the known results for Q^f {y) and pj)^'^ {y)l It turns out that two communities 
have derived different formulae for these very same quantities using different languages: 
HFMA on the one side, and determinants or Pfaffians of Bessel functions on the other 
side (see two next subsections). It is one of the goals of this paper to highlight this 
connection, which is not widely appreciated, and to actually prove the equivalence of 
these formulae by simple algebraic methods. As a result of this (and of the previous 
section) we will derive some new compact results in the second language in cases where 
only HFMA or no results were previously known. 

We should also mention here that a third way exists to compute the cumulative 
distribution. Formally it is given by the Fredholm determinant of the corresponding 
Bessel kernel. This has been mapped to expressions containing a transcendent of 
the solution of the Painleve V for /3 = 2 [30] and /? = 1,4 [H]. In particular, the 
distributions for different /3s can be related to each other. However, this relation 
does not allow us to generate more explicit expressions in terms of determinants and 
Pfaffians below. 



4.-1- Hypergeometric functions of matrix argument (HFMA) 

In order to establish the claimed equivalence we first need to state the results in the 
first language. Forrester [55] and Chen et al. [53] independently gave expressions for 
P!)^\y) and Qi^\y) valid for ah /3 > and integer m = (/3/2)(i/ + 1) - 1 e N, in terms 
of HFMA as: 

qW) (y) = e-/3y/8 ^piPm (_. 2m/ (3; {y/A)l„,) , (4.8) 

P(«(y) = A™,^ y"e-^^/8 qF^^/'^-; 2 + 2m//3; (y/4)I™) , (4.9) 
where the constant A^.p is given bjjjj 

and Im is the identity matrix m x m. For definitions and properties of HFMA we refer 
to Appendix |A] We note that for = 1 integer m implies that v — M — N must be 
odd. 

The formulas above were specialized in [23] to the two cases (a) m = 1 for all /3 > 0, 
and (b) for /? = = 2, yielding respectively: 

(a) Q:,l\4/,)-i(2/)-2-i+^/''r(2//3)e^^^/«yV2-i//.I,_^(^) (4.n) 

(b) Q[^=2\y) =e-y''{H^f~M^f) (4.12) 

* We note in passing that the large y S> 1 behavior of Qi,^'(j/) is known from |40l . see also |42l for 
the calculation of leading term using a Coulomb gas method. 

fj Note that the expression for the constant A^^p, eq. (2.15c) in 1281 (and subsequently in [23]), as 
well as eq. (|4.11| as it appears in [23) contain misprints. 
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where lp{x) is a modified Bessel function of the first kind. The appearance of Bessel 
functions for special instances of formulas ( 4 
will see now. 



and (4.9 1 is not at all accidental, as we 



4-2. Determinants and Pfaffians of Bessel functions. 

In the context of Random Matrix Theory applied to effective theories of Quantum 
Chromodynamics (QCD) or QCD-like theories, the large- distribution of the smallest 
Dirac operator eigenvalue has been studied intensely for the symmetry classes /3 = 1,2 
(and 4), see [TU] for references. It turns out that the Dirac eigenvalues (occurring 
in ± pairs due to chiral symmetry) are precisely distributed according to the WL jpdf 
(2.8) upon the mapping yf — Xi. The relation between the smallest Dirac eigenvalue 



^pl?^(s), and the smallest WL eigenvalue distribution trivially follows: 

^(«(.) = 2.P(«(.^) 
0(«(.) = Q(«(.2), 



where we also give the relation to 
Dirac picture. 



We still have *pi^^(s) 



the cumulative distribution Cl[f\s) 



(4.13) 
(4.14) 

in the 



~£li^'^'{s) and the standard normalization 



ds (s) = 1 is ensured. 
The following results hold for our symmetry classes: 



fn(/3=l) 



2 



s)e 
^s)h{s) 



/«Pf[(z-J-)W3(s)], 
/8-s/2 



-m+i 



m e N, 



Ms) 



(4.15) 
(4.16) 
(4.17) 
(4.18) 



where Pf(yi) ~ ^/detA stands for the Pfaffian of the antisymmetric matrix A. 

Most of these results were known in the literature, apart from the normalization 
constant Cm given in eq. (4.41), and the new eq. ( 4.18[) to be derived below. For 
(3 = 2 the distribution eq. (4.15) was derived in [26l 1341 135]. For /3 = 1 and ly — eq. 
(4.17) was shown in [57], and then extended to all odd in up to normalization. 

We note in passing that for /3 = 4 the distribution for v = is known explicitly from 
[?7] (see eq. (4.11) and take its derivative), but the results quoted in [25| for /3 = 4 
only hold for ni — I3{v + l)/2 — 1 even, that is for half-integer values of v. Finally we 
mention that for small v all distributions for (3 = 1,2,4 can be very well approximated 
by using finite N = 2 results [13], like in the Wigner surmise. 

Turning to the corresponding cumulative distributions the following results hold: 

0(^=2)(s)= e-^^/4det[I,_,(.)]^_^.^^ ..^^ , ven (4.19) 
2 [35j , and for /3 = 1 and odd v — 2m + 1 our new result reads 

.=2™+i(5) = 2"s-"e-^'/«Pf [(j - k)h+k+,is\k=^m+^^,...^,n-i € N. (4.20) 

Before providing a general algebraic proof of the equivalence of the two languages in 
subsections |4.4| and |4.5[ we can quickly check the agreement between the two different 
formulations above, starting from the special HFMA cases in eqs 
eq 



for /3 



Using the map (4.14) 



in eqs. (4.19) and (4.20) respectively. 



(4.19) obviously contains eq. 



Specifying /? in eq. (4.11) leads to the special cases (/? = 2, 



(4.11) and (4.12) 



(4.12) as a special case 
' = 1) and i/3 = l,v 



3) 



Additionally one can differentiate these and arrive at corresponding special cases in 



eqs. (4.15) and (4.17) respectively, after using some Bessel identities. 
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4-3. Results for the scaled smallest eigenvalue distribution for P — 1 and v = 0, 2. 



In this subsection, we give two results for the scaled smallest eigenvalue distribution 
at even v — Q and v = 2. In principle the universality eq. (4.4), that is the agreement 



between the WL and FTWL distributions in the microscopic large-iV limit, has not 
been shown for fi — \ with even although there is little doubt that it extends from 
the odd V case. We explicitly check here that this is indeed the case by taking the 
microscopic large- iV limit for FTWL at = starting from the known eq. (3.201 for 
finite- TV and getting to the same result as in WL. For v — 2^e derive a new result that 
was unknown even for the microscopic limit of the unconstrained WL ensemble. 

For 1/ = we first recall how the known microscopic result for WL is obtained. 
Starting from the finite- iV expression eq. (3.9) together with eq. (3.14) we need to 



know the asymptotic limit of the Tricomi confluent hypergeometric function, in the 



scaling limit eq. (4.2). It was derived in Corollary 3.1 of (Ml, alternatively it follows 
from eq. 13.3.3 [S]: 



lim 



V2 



1 U 



X 

2n 



Collecting all prefactors we thus obtain 



1 



2 



(4.21) 



(4.22) 



The relation eq. (4.131 maps this to the eq. (4.17). This result has to be compared 
with the limit eq. (4.1 1 of cq. (3.20) in FTWL ensemble. Because the limit of the 



hypergeometric function is more involved we defer the derivation to AppendixjBj finding 
complete agreement as in eq. (4.4 1: 

1 



j(/3=l) 



(y) 



u=Q (y) = lim 

5-OC 



1 

4/V3 



(FT 

Pn 



V 
4Ar3 



(WL) 



(y) = lim 



AN 



(WL) / J^A 



Pn,u 



(4.23) 



This extends the expected universality to the case of even v 
For ly 



0. 



2 the corresponding microscopic limit of (3.21) for FTWL is already rather 
cumbersome, involving the asymptotic of a sum of hypergeometric functions. We 
therefore restrict ourselves to compute the corresponding microscopic limit in the WL 
ensemble - which is also new - and conjecture that the universal relation eq. (4.23) 



extends to v = 2 and in fact to all higher even v^s. Combining eqs. (|3.9| and (|3.15) 
we have to compute 

' N 



AN . 



-(2) 



U 



lim 

1 

'^2' 



= ^/2^ \an) 

8N J 



'V/S 



-y \ N-_ 

ml 4 



-u 



1 



8N / 



Eq. (4.21) together with the Laguerre asymptotic for negative argument, 
lim m-^Ll^^-x/m) = x'''^\{2yG) , 

m— >-C30 

yields the following final answer, 



(4.24) 



Using eq. (4.13) this is mapped to eq. (4.18) as claimed. 
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4.4- Equivalence proofs for (3 = 2 

In this subsection and in the following, we provide an algebraic link between the HFMA 
and the Bessel determinant languages in the spirit of earlier works [HJ |3S] . We start 
from the formula in eq. (4.9 ) which expresses the scaled smallest eigenvalue distribution 
(with or without fixed trace constraint) in terms of HFMA: 



(4.26) 



Next, we use the following integral representation due to Forrester [2F valid for integer 
values of A and c: 

(c-l)m/2 



oi^f/^' (-;c+^(,7i-l);rf„ 



where we have defined 



Bm{c, A) 



1 

2^ 

cos 6j i{c—l)9j 



r(l + A/2)r(c + A(j-l)/2) 

r(i + Aj72) 



A„(c,A) = [] 

and the Vandermonde determinant of the angles 
A(e'^^) = ]^(e'^'= ■ 



A(e"') , (4.27) 



(4.28) 



(4.29) 



Comparison between (|4.27|) and (|4.26|) yields c = 3 and A = 2. Denoting by * complex 

(4.30) 



conjugation, we can use the Andreief identity |46j 

dOe'^^'-'^^Loie) 



]J(i0ja;(6lj)A(e"')A*(e'®) =m! det 

j — \ L-^ J £,A: — l,...,m 

We thus obtain for the m-fold integral in the second Une of eq. ( 4.27[ ) 

m! det [ r d0e'«(^-fe+2)e2yscosel ^ ^, ^^^y. [l,_fe+2(2V^)]^^^^^ 

(4.31) 

where we have used the following integral representation for the modified Bessel 
function (valid for integer index only): 

1 



dt 



' cos(ni) , n e N 



(4.32) 



Setting X = y/A in (4.31) and simplifying all prefactors, we eventually get from (4.26): 

P^^'Jrniy) = le-y/^det [Iw+2(W)].,,=i...... , ^ e N . (4.33) 

This is identical to eq. ( |4.15 1 after switching to the Dirac picture eq. (4.13). 

The proof for the cumulative distribution goes along the same lines, merely changing 
the coefficient in front and the index shift of the Bessel function. We obtain the 
following from eq. (4.8) by identifying c = 1 and A = 2 in the representation eq. (4.27) 
and collecting all prefactors: 

Q^f=^Hy) ^ e-y^\Fi'\-;m;{y/4)lr^) = e'y/^ det[h_,{^)l^^^^_^ ,z.eN(4.34) 

This corresponds to eq. (4.19) using the map eq. (4.14). Note that it is highly non- 
trivial to derive eq. (4.34) from eq. (4.33) using eq. (4.5). 
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4-5. Equivalence proof and new results for l3 — 1 



We start again from the formula in eq. (4.9) expressing the scaled smallest eigenvalue 
distribution (with or without fixed trace constraint) in terms of HFMA: 

P!.itl+i{y) = Ans y'^e-y/^ oF^'/^^ (-; 2m + 2; , m e N . (4.35) 



Next, we use again the integral representation in eq. (4.27) and x = y/A to write: 



(1/2) 



(-;2 + 2to;xI„0 =B„(4,4) 



3m/2 



(4.36) 



Comparison between (|4.27|) and (|4.35|) indeed yields c = 4 and A = 4. Next, we use 

(4.37) 



section 11.5 in [5], combining eqs. (11.5.2) and (11.5.4) there to write 



j<k 



= det [e'P^^" pe'P^^"] 



where the determinant on the right hand side is of size 2m x 2m and indices run as 
follows, 1 < j < m and p = -(m- 1/2), -(m - 3/2), . . . , (m - 3/2), m- 1/2. 

Now we can define the two sets of functions {4>k{0)} — e'*^^ and {V'fe(^)} — ke^^^ and 
thus use the de Bruijn identity [47 : 



<j<m 



/lit 
\Y dXjUj{Xj) det [(t)i{Xj) Mxj)]i<,<2m.i< 

dx Lo{x) {(l)i{x)^j{x) - (l)j{x)tpi{x)) 



m! Pf 



(4.38) 



l<'ij<2m 



to evaluate the m-fold integral in the second line of eq. (|4.36|): 
ml Pf 



de e2yScoseg3ie(gifce^.gye _ ^ijOf^^ike^ 

= m!(27r)'" Pf [{j - k)h+k+j i^V^)] (4.39) 

where the indices {k,j) all run from — (m— 1/2), —(m — 3/2), . . . , (m ~ 3/2), (m — 1/2), 
and we have used again the integral representation for the Bessel function (4.32). 
Setting X = j//4 in (4.39) and simplifying all prefactors, we eventually get for eq. 

(4.40) 



(4.35): 



^,£2i^i(y) = C„,y-"/2e"«/«Pf [(j- - fc)l3+fe+,(V^)] 
where is the normalization constant: 
y^(2m + 1)!! 



C'ri 



16 r(m + 3/2) 



(4.41) 



Applying now eq (4.13), we get the distribution in the Dirac picture (4.16) including 
its explicit normalization constant for any v, which was previously unavailable. 

The proof for the cumulative distribution easily follows from (4.8), and we only 
quote the result which is new, after having identified c = 2 and A = 4 in eq. (4.36) 
instead: 

Qi^=2l+iiy) = ^-"'^ oF^'/'H-, 2m; (jy/4)I„) 

= 2"^y'"^^'e-y^'Pf[ij-k)h+k+,iVy)] ■ (4-42) 
The range of indices is the same as in eq. (4.39). Once again, deriving eq. (4.42) 
directly from eq. (4.40) using eq. (4.5) is highly non-trivial. 
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4-6. Moments for 



N 



We point out a universal expression for the large N decay of the £th moment of the 
smallest eigenvalue, valid for both (3 = 1 and 2: 



for — J- oo 



(4.43) 



where we have defined the universal coefficient 

/"OO 

JO 

The scaling with different powers of N trivially follows from the different spacings 



(4.44) 



in the definition of the microscopic limit eqs. (4.1 1 and (4.2), as we will show now 



Starting from the definition of the average of the smallest eigenvalue for FTWL (2.18) 
we have 

'^^'^-"'(^^=(4wi 4M^- (4I3) (4-45) 



after making a change of variable x — y/AN^. This implies, using the limit (4.1 ) 

"(4A^^)^(ALn)Ki 



lim 



dy /Pi^\y) 



(4.46) 



The same argument for WL with a different change of variables, x = y/4:N, leads to 



lim 

TV— foo 



(4A^)^(AL„)K 



dy /P^^Hy) ■ 



(4.47) 



Alternatively the universal right hand side can be expressed in the Dirac picture, using 



the map eq. (4.13), as it is given in eq. (4.44) 



100 120 



Figure 4. The behavior of {X^ij^)^^'^\ for increasing A' and for u = 1,3 (left 
and right panel respectively). The curves converge for Af — 00 to the limiting 
values (1/4) = 2 (left) and (1/4) ki,3,i = - 1 6.38906 (right). 



Let US I 

integral 



jive a few examples, by simply inserting eqs. (4.17) and (4.15) into the 
1 

= £ / dss2^(2 + s)e-^'/8-^/2 

= 2^^ \-y/2 IT {1/2 + I) iFi (1/2 + £; 3/2; 1/2) 

+r (1 + I) (iFi (1 + t- 1/2; 1/2) - iFi (1 + t- 3/2; 1/2) )] (4.48) 

1 

= - ds s^'+^e''"'^ = 4^r(l + £) . (4.49) 
2 Jo 
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Specializing these results to the case of the first moment £ — I (average value), 
obtain: 



we 



K£=l,,y=0,/3 = 1 



K£=l,i/=0,/3=2 



0.688641. 



1 



(4.50) 
(4.51) 



in complete agreement with |1H I18j (after switching to their conventions). However, 
the general asymptotic relation ( |4.43 1 provided here allows to derive analogous, new 
results for any v for which the smallest eigenvalue in the Dirac picture is known. 

We have numerically verified the scaling behavior above for FTWL ensemble for the 
two cases = 1,3 in fig. |4j We plot the combination A^'^(Amin)jv,i/,i as a function of 
A'^ for = 1, 3 (left and right panel respectively), which nicely converge for TV — )■ oo to 



the predicted values (1/4)ki^i_i 
used (|4.44|) along with (|4.16). 



2 and (1/4)ki 



3,1 



I « 6.38906, where we have 



5. Conclusions 



In this paper we addressed the smallest eigenvalue distribution in fixed-trace Wishart- 
Laguerre (FTWL) ensembles of random matrices, as well as its cumulative distribution 
and moments, and the corresponding quantity for unconstrained Wishart-Laguerre 
(WL) ensembles. Our motivation comes from the statistical description of entangled 
random pure states in bipartite systems of size MN , where the distribution of 
the smallest Schmidt eigenvalue provides useful information about the degree of 
entanglement of these states, and as such it has been subject to intense scrutiny in 
recent years. 

In the first part, we derived new results for the FTWL ensemble with orthogonal 
symmetry. Here the constraint plays an important role, and we computed explicit 
expressions for the smallest eigenvalue by Laplace inverting a recursion relation given 
by Edelman for the unconstrained WL ensemble. In particular our results extend very 
recent expressions to even values oi v = M — N . Examples were given for several even 
and odd values of v and checked via numerical simulations. 

In the second part, we studied the microscopic large- TV limit at fixed both for 
orthogonal and unitary ensembles. Here the constraint (after a proper rescaling) is 
immaterial due to universality. We proved the equivalence of two different sets of 
results in the literature, given in terms of either hypergeometric functions of matrix 
argument (HFMA), or in terms of Pfafhans and determinants of Bessel functions. Our 
technical tool was to translate HFMA into their multiple integral representations and 
to use the classical Andreief-de Brujin identities. Apart from this equivalence proof, 
our new results include the cumulative distributions, the universality of the smallest 
eigenvalue for = and an explicit expression for = 2, all in the orthogonal symmetry 
class. In addition we have determined the precise asymptotic behavior of all moments 
for large TSS ^ in both symmetry classes. 

Several open problems persist, in particular how to find a closed universal formula 
for all even values of v in the orthogonal symmetry class in the microscopic limit. 
Furthermore, it would be interesting to extend these results to the symplectic symmetry 
class. Here a compact expression is only known for the special integer = 0. The 
problem is the apparent lack of a suitable integral representation for HFMA in order 
to proceed. 
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A. Hypergeometric functions of matrix argument (HFMA) 

The hypergeometric function of a matrix argument takes a complex symmetric matrix 
(m X m) X as input and provides a real number as output. More precisely, let p > 
and q > he integers, and {xi, . . . , Xm} be the real eigenvalues of X. We have: 

(ai, . . . , a,; 6i, . . . , 6,; X) ^ ^ luu.m lm ^''^ ^^^^ 

where > is a parameter, the symbol n \- k means that k = (ki,K2,...) is a 
partition of k (i.e. ki > ^2 > . . . > are integers such that \n\ — ki + ^2 + . . . ~ 

fc), and the following symbol (a)if'' — 11(1 j)GK ('^ ~ ^ + i ~ ^ generahzed 



Pochhammer symbol. In (Al), cli\X) is a Jack polynomial, i.e. a symmetric. 



homogeneous polynomial of degree |k| in the eigenvalues {xi} of X [48l |4^. Jack 
polynomials generalize the Schur function, the zonal polynomial and the quaternion 



zonal polynomial to which they reduce for j3 = 1,2,4 respectively. The HFMA (All 
generalizes the ordinary hypergeometric function to which it reduces when m — I. 



The series (Al) converges for any X if p < g; it converges if max^ \xi\ < 1 and 
p = q + 1; and diverges if p > g + 1, unless it terminates. An efficient evaluation of 
HFMA is now made possible by an algorithm devised by Koev and Edelman [50], which 
we used extensively to check numerically all the results given in terms of HFMA in this 
paper. 

B. Microscopic limit of p^^Lo(^) ^'^^ /3 = 1 



In this appendix we take the microscopic limit eq. (4.1) of the smallest eigenvalue for 



hm - . ( ^ 1 = hm NnNnNy2) ( y ^ 



1 and in FTWL eq. ( |3.20D 

J_ (FT) ( JJ_\ _ 

N-^^oo 4iV3Pw.^=o Ua^^; j^-;^ Ar32W+ir(iV/2)r((Ar2 +N- 2)/2) \AN^ 

X 1 %\ iFxX — ; \N . Bl 

While the first factor in the second line yields an exponential prefactor, 

hm (l-^y = e'y^\ (B2) 

the asymptotic limit of the hypergeometric function requires more care. Using eq. 
9.132.1 in [21] we have that 

/A^ + 2 N-1 + 
^ ^ V 2 ' 2 ' 2 ' y 



_ N + 2 



r^F + f-i r(-i) 



X^^^^ At, ^.Fi ( ^ + 1, ; ■,-;X 



rff-i)r(^-2K 'V2 ' 2 '2 



A^ -1 5 



_ N + 1 



r T + r(i) 1 



2 2 

with 



^ ' r(f + i)r(-^' — rv2^i (^^' — "2;"9;^) (^^^ 



X = [ — - N + 1] . (B4) 
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In order to apply the standard asymptotic limits, eqs. 9.121.9-10 in 

: cosh(z) 



/ 1 z2 
J^oo^^i(,^'^'2'4^ 



sinh(z) 



(B5) 
(B6) 



where we identify j = N/2 and k = N'^/2, we still have to shift the third index of the 
two hypergeometric functions in eq. (B3|. Using eq. 9.137.1 in [51j we obtain for the 
first 



■ N 

, 2 2 



N 

y 
-1 



' 2 ' 2' 



-2)X 



3 

\ 

N 3 
~2 ^ 2 



[X - 1) 2F1 



N 

y 

N 

y 



' 2 ' 2' 



X 



(B7) 



2^1 



and for the second 
'A^-1 

y^' ^ 
-1 



2;~-:X 
' 2' 



N 

1 

2 



2 I ^ 2F1 



N-1 
2 ' ^ 



2;^;X 
' 2' 



AT 

^ + y 



^ 2^1 



iV- 1 
2 ' ~Y 



2;l;X 



(B8) 



Collecting all the results from above and all factors we obtain as a final result 



cosh 



2 



: sinh 



— 1 cosh 1 I — sinh 



2 



2 



(B9) 



which agrees with eq. (4.22) as we have claimed. 
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